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Abstract 

The long-time asymptotics is analyzed for finite energy solutions of the ID discrete 
Schrodinger equation coupled to a nonlinear oscillator. The coupled system is invariant 
with respect to the phase rotation group U{1). For initial states close to a solitary wave, 
the solution converges to a sum of another solitary wave and dispersive wave which is a 
solution to the free Schrodinger equation. The proofs use the strategy of Buslaev-Perelman 
[2]: the linerization of the dynamics on the solitary manifold, the symplectic orthogonal 
projection, method of majorants, etc. 
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1 Introduction 



Our main goal is the study of the distinguished dynamical role of the "quantum stationary 
states" for a model f/(l)-invariant nonlinear discrete Schrodinger equation 

i7p{x,t) = -AL'^{x,t) -6{x)F{ilj{0,t)), xeZ. (1.1) 

Here F is a continuous function, 6{x) = 6qx and stands for the difference Laplacian in Z, 
defined by 

AliIj{x) = ipix + 1) -2iIj{x) +iIj{x -1), xgZ 

for functions ip : Z ^ C Physically, equation fll.ip describes the system of the free discrete 
Schrodinger equation coupled to an oscillator attached at the point x = 0: F is a nonlinear 
"oscillator force" . 

We identify a complex number = i/ji + i'il>2 with the real two-dimensional vector \l/ = 
(V'1,'02) e and assume that the vector version F of the oscillator force F admits a real- 
valued potential, 

F(^) = -Vt/(^), *eM^ f/GC^(M^). (1.2) 
Then (II. ip is a Hamiltonian system with Hamiltonian 

U{m) = ^(-A,.^, + f/(^(0)) = ^(Vl^, Vl^) + f/(^(0)), (1.3) 

where (■) stands for the scalar product in /^(Z), and Vl'?/'(x) = tp{x + 1) — ip{x). 
We assume that U{ip) = u{\ip\'^) with u G C^{R). Therefore, by 

F(V') = adVHV', V-eC, aeC\R), (1.4) 

where ad^/'H is real. Then F(e^V) = e'^Fitp), 6 G [0, 27r] and F(0) = for continuous F. 
Hence, e^^ip{x,t) is a solution to (11. ip if ip{x,t) is. Therefore, equation (11. ip is f/(l)-invariant 
in the sense of [5], and the Nother theorem implies the conservation of norm: 

\\m\\ = \\mi 

Here and below || ■ || = || ■ 

The main subject of this paper is an analysis of the special role played by solitary waves or 
nonlinear eigenfunctions, which are finite energy solutions of the form 

ip{x,t) =i;^{x)e''^\ cuER. (1.5) 

The frequency u and the amplitude ipuj{x) solve the following nonlinear eigenvalue problem: 

- = -Al^^(x) - 6{x)F{^^{0)), xeZ (1.6) 

which follows directly from (II. ip and (11.40 since G M. The solitary waves constitute a two- 
dimensional solitary manifold in the Hilbert phase space /^(Z). 
We prove the asymptotics of type 

^i; t) ~ ^.±e^"±* + (t)$±, t ^ ±oo, (1.7) 

where W(t) is the dynamical group of the free Schrodinger equation, ^± G /^(Z) are the 
corresponding asymptotic scattering states, and the remainder converges to zero as (9(|t|~^/^) 
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in global norm of /^(Z). The asymptotics hold for the solutions with initial states close to the 
stable part of the solitary manifold, extending the results of [H [21 [Sj [101 [HI [12] to the equation 

For the first time, the asymptotics of type fll.Tp were established by Soffer and Weinstein 
[T^ [T5] (see also [13]) for nonlinear f/(l)-invariant continuous Schrodinger equation with small 
initial states if the nonlinear coupling constant is sufficiently small. The next result was obtained 
by Buslaev and Perelman |2] who proved that the solitary manifold attracts finite energy 
solutions of a ID nonlinear f/(l)-invariant translation invariant Schrodinger equation with initial 
states sufficiently close to the stable part of the solitary manifold. For a more lengthy discussion 
of our motivation, and of previous results in the literature ([3l HH [121 [131 El HSl [IE]) we refer 
the reader to the introduction of [Ij. The asymptotics of type (11.71) for nonlinear discrete 
Schrodinger equation are obtained for the first time in the present paper. 

Let us note that we impose conditions which are more general than the standard ones in the 
following respects: we do not hypothesize any spectral properties of the linearized equation, 
and do not require any smallness condition on the initial state (only closeness to the solitary 
manifold). This progress is possible on account of the simplicity of our model which allows an 
exact analysis of spectral properties of the linearization. 

Let us comment on the general strategy of our proofs. We develop the approach [H |2l [6l [8] 
for our problem. Firstly, we apply the symplectic projection onto the solitary manifold to sep- 
arate the motion along the solitary manifold and in transversal direction. Secondly, we derive 
the modulation equations for the parameters of the symplectic projection, and linearize the 
transversal dynamics at the projection of the trajectory. The linearized equation is nonau- 
tonomous, and this is one of the fundamental difficulties in the proof. This difficulty is handled 
by the introduction of an autonomous equation (by freezing the time) with an application of 
the modulation equations to estimate the resulting additional error terms. A principal role in 
the rest of the proof is played by the uniform decay of the frozen linearized dynamics projected 
onto the continuous spectrum, and the method of majorants. 

The paper is organized as follows. In ^we describe all nonzero solitary waves and formulate 
the main theorem. In ^ we enumerate some properties of the linearized equation. In §1] 
we establish the time decay for the linearized equation in the continuous spectrum. In [|5] 
the modulation equations for the parameters of the soliton are displayed. The decay of the 
transverse component is proved in §21 and §71 In §H]we obtain the soliton asymptotics (II. 7p . In 
Appendix we study the resolvent of linearized equation. 

2 Solitary waves and statement of the main theorem 

Our main results describe the large time behavior of the global solutions whose existence is 
guaranteed by the following theorem 

Theorem 2.1. J3 Theorem 3.1] 

i) Let conditions / (1.2|) and MA\) hold. Then for any i/jq & l"^ = l'^{'L) there exist a unique 
solution ip G C;,(M, /^) to the equation ( fl.lp with initial condition ip{x,0) = ipoi^)- 
a) The value of energy functional and the norm of solution is conserved: 

W(t))=Wo), ||^(t)|| = ll^oll, tER. (2.1) 
In [9] we give a complete analysis of the set of all solitary waves. There exist two different 
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sets of nonzero solitary waves. The first set 5+ corresponds to G (0, oo) and the second set 
S- corresponds to G (— oo, —4). 

Denote by k{u) the positive solution of the equation cosh/c = \u! + 2|/2. 

Lemma 2.2. /P| Lemma 4-^]) The sets of all nonzero solitary waves is given by 

S+ = |^/'^e*^ = Ce^^-^(")l^l : u G (0,cx)), C > 0, sinhA;(cu) = a{C^)/2 > 0, 6 e [0,27r]} 

5_ = |^/>^e*^=C(-l)l^le^^-^(^)l^l: G (-oo, -4), C>0, sinh A;(cu) = -a(C2)/2 >0, eG[0,27r]} 

Corollary 2.3. T/ie setS+ resp. S- is a smooth manifold with the co-ordinates 6 mod 27r 
and C >Q such that aiC"^) > resp. aiC"^) < 0. 

In the case of polynomial F the condition on C means that C is restricted to lie in a set 
which is a finite union of one- dimensional intervals. The value C = corresponds to the zero 
function iIj^^^x) = which is always a solitary wave since F{0) = 0, and for uj G [—4, 0] only the 
zero solitary wave exists. 

In ([9]) we proved that the parameters 9, oj locally also are smooth coordinates on S± at the 
points with a'(C^) 7^ 0. We will also need the following result from [9]. 

Lemma 2.4. For C > 0, a > 0, and a' 7^ (4a + a^)/(4C^) we have 

[ \M^)\^dx^O. (2.2) 



The soliton solution is a trajectory 4'oj{t)ix)e^^^^\ where the parameters satisfy the equation 
6 = uj, uj = 0. The solitary waves e^^il)^{x) map out in time an orbit of the U{1) symmetry 
group. This group acts on the phase space /^(Z) preserving the Hamiltonian Ti. 

The real form of the solitary wave is e-'^^^ where = {1p^^{x), 0), and j is the 2x2 matrix 

Linearization at the solitary wave e^^^^ leads to the operator (cf. [T, "3]) 

Di C 
D2 

where Pi is the projector in acting as 



B = -AL + u-SixMC') + 2a'{C')C'Pi]= ( ^ ), (2.4) 



X2) [0 

Di = -Al + UJ - 6{x)[a + 2aC\ D2 = -A^ + uj - 6{x)a. 

Let C = j^^B. We will show in Appendix A that the continuous spectrum of C coincides with 
C_ U C+ = + 4), —ioj] U [ioj, i{uj + 4)]. The point belongs to the discrete spectrum, and 

the dimension of its invariant subspace is at least 2. If a' 7^ (4a + a^)/(2C^) then the invariant 
subspace associated to the eigenvalue A = is of dimension exactly 2. We will analyze only 
the solitary waves with a' ^ {0; (4a + a'^)/(2C^)}. We assume also more specific condition 
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Definition 2.5. We say the solitary wave ipi^{x)e^^ , satisfies the spectral condition {SP} if 

1) a'^{0;(4a + a3)/(2C2)} 

2) there is no eigenvalue except A = 0. 

In Appendix A we give an example wlien tlie condition {SP} liolds. If condition {5'P} is 
true for a fixed value cuq, it is also true for values in a small interval centered at ujq. The 
condition {SP} ensures orbital stability of solitary waves. 

The functional spaces we are going to consider are the weighted Banach spaces = /^(Z), 
p G [1, cxd), /5 G M of complex valued functions with the norm 

\\u\\f^ = \\{l + \x\fu{x)\\i.. (2.5) 

Let us denote by W{t) the dynamical group of the free Schrodinger equation. Our main theorem 
is the following: 

Theorem 2.6. Let conditions / II. 2)) . and ( (1.4)) hold, (3 > 2 andilj{x,t) G C(M.,P) be the solution 
to the equation ( (l.lj) with initial value ipoix) = ip{x, 0) G /^fl/^ which is close to a solitary wave 

d:=Uo-iJ^,e''%2r,ii<^l. (2.6) 

Assume further that the spectral condition {SP} holds for the solitary wave with u = ujq. Then 
for d sufficiently small the solution admits the following asymptotics: 

t) = ^^J^^' + Wif)^^ + r±(t), t ^ ±oo, (2.7) 

where $± G Z^(Z) are the corresponding asymptotic scattering states, and 

||r±(t)|| = 0(|t|-V2), t-.±oo. (2.8) 



3 Linearized evolution 

In this section we summarize the properties of the linearized evolution which will be needed. 
The proof of these properties one can find in Appendix A and in [Tj. The linearized equation 

reads 

x(a;,t) = Cx(x,t), C := r'B = f ° V (3.1) 



Theorem 12.11 generalizes to the equation fl3.ip : the equation admits unique solution x{x,t) G 
Cb(]R, /^) for every initial function x(^) 0) = Xo ^ Denote by e^* the dynamical group of 
equation (13. ip acting in the space Then (12. ip implies that 

lle'^^Xoll = IIXoll, teK. (3.2) 

The resolvent R(A) := (C — A)^^ is an integral operator with matrix valued integral kernel (see 
Appendix A) 

R-(A, x, y) = r(A, X, y) + P(A, x, y), (3.3) 



where 



T{\,x,y) 



^\ 

4 sin 01 4sin0_ 



I 



\ 4 sin 6'. 4sin0_ / 



' ^ie+\x-y\ _ ^i9+{\x\ + \y\) _l(^QiS+\x-y\ _ ^ie+{\x\ + \y\)Y 



^ie.\x-y\ _ ^ie_(\x\ + \y\) ^f^i0.\x-y\ _ ^ie.(\x\ + \y\)\ 



(3.4) 



Here 9±{\) is the root of 2cos6'-|- = u + 2 ± iX defined with cuts in the complex A plane so 
that 9±{X) is analytic on C\C±, and Im6'±(A) > for A G C \ C±. The constants a, j3 and the 
determinant D = D{X) are given by the formulas 

a = a + a'C'^, (3 = a'C'^, D = 2ia{sm 0+ + sin 6'_) - 4 sin 6*+ sin 6'_ + - 

The poles of the resolvent correspond to the roots of the determinant D{X). 

If spectral condition {SP} holds then the determinant has the only root A = with the 

multiplicity 2. 

Observe that fll.6p and its derivative in uj give the following identities: 

B2i^^ = Di(9^^^) = (3.6) 

These formulae imply that the vectors j^^j and d^^^^j lie in the generalized two dimensional 
null space of the non-self-adjoint operator C defined in (13. ip and 

Cj$^ = Cd^<^^ = (3.7) 

The symplectic form Q for the real vectors ip and rj is defined by 

fi(V',r/) = {ifjjT]) = ^{ipir]2-'^2Vi)- (3-8) 

By Lemma [2.41 

^^(j$., d^^J) = J \ipj^dx ^ 0. (3.9) 

Hence, the symplectic form Q is nondegenerate on X°, i.e. is a symplectic subspace. 
Therefore, there exists a symplectic projection operator P° from onto X° represented by the 
formula 

= /eh Icfy A i^^jd.'^^h^^ + {^,'^^)dA] (3.10) 

Denote by = 1 — P° the symplectic projector onto the continuous spectral subspace. 

Remark 3.1. On the generalized null space itself C"^ = by ( 13. 7j) . and so the semigroup e**" 
reduces to 1 + Ct as usual for the exponential of the nilpotent part of an operator. 

4 Time decay in continuous spectrum 

Due to Remark [3.11 we see that the solutions xit) = ^^^Xo of the linearized equation (13. ip do 
not decay as t — oo if P°Xo 7^ 0. On the other hand, we do expect time decay of P^x(^)) as a 
consequence of the Laplace representation for P'^e^^: 



P^e 



c^Ct 



— — [ e^*(R(A + 0) -R(A-O)) c/A. (4.1) 
27ii J ^ 



2m 

c+uc 



The decay for the oscillatory integral is obtained from the analytic properties of R(A) for 
A G C+ UC_. 
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In Appendix A we will show that -D(A) 7^ for A G C+ U C_. Clearly in order to understand 
the decay of P'^e**", it is crucial to study the behavior of R(A,x, ?/) near the branch points 
A = ±iuj and X = ±i{uj + A) (where sin^± vanish). 

We deduce time decay for the group P'^e**^ by means of the following version of Lemma 10.2 
from [7j, which is itself based on Zygmund's lemma [IHI p. 45]. 

Let : [a,b] — >■ B be a function with values in a Banach space B. Let us define the 
B- valued function 



I{t) = J e-'^'J^iu) du. 



Lemma 4.1. Suppose that T{a) = T{h) = 0, and T" G L^(a + 5, 6 — 5; B) /or some 5 > 0. 
Moreover, 

in the norm of B for some p G (0, 1). Then I{t) G Cb{a + £, 6 — £; B) for any £ > 0, and 

I{t) = 0{t-^'P) as t^ 00 in the norm of B. 

We will apply Lemma 14.11 to the function J^(A) = R(A + 0) — R(A — 0) with values in the 
Banach space B = B{lj^, Z^^) , the space of continuous linear maps Ij^ for any {3 >2. 

Theorem 4.2. Assume that the spectral condition {SP} holds so that \ = Q is the only point 
in the discrete spectrum of the operator C = C{uj). Then for (3 >2 

||p-e^*||g = 0(^3/2), t^oo. (4.2) 
First we use the formulas (14. ip and (13. 3p to obtain 

-27riP^e<^*= j e^*(r(A + 0) -r(A-0))tiA + j e^*(P(A + 0) - P(A - 0)) dA (4.3) 

c+uC- c+uc_ 

Next we apply Lemma 14.11 to each summand in the RHS of (14. 3 p separately. Then Theorem 
14.21 immediately follows from the two lemmas below. 

Lemma 4.3. If the assumption of Theorem \A.2\ hold then 

j e^*(r(A + o)-r(A-o)) rfA = o(r3/2^, ^44^, 

c+uC- 

in the norm B. 

Proof We consider only the integral over C+ since the integral over C_ can be handled in 
the same way. The points \ = iu and A = z(u; + 4) are the branch points for 9^, therefore, if 
A G C+ then since 0- is continuous across C+ 

r(A + 0) - r(A - 0) = r+(A + o) - r+(A - o), 
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where F"*" is the sum of those terms in F which involve 9+. Let us consider, for example, F^^. 
The expression (13 ■4p implies for ?/ > that 



f 0, x < 



rii(A,x,y) = < 



4 sin 9+ 

^ie+xf^-ie+y _ ^ie+y\ 



< X < y, 



X > y. 



4 sin 9+ 

For A G C+, the root 9+ is real, and 9+{X + 0) = -9+{\ - 0). Then, for y > 0, 

sin 6'+|x| sin6'+|y| 



F+(A + 0,a;,i/)-F+(A-0,a;,y) = -9(0;)- 
where 0(x) = 1 for x > and zero otherwise. Recall that 



sin 6* J 



(4.5) 



6*+ = 9^{i') = arccos — , sin 6'+ = - A/(^^Z7+Ty(T+~c7^^^, u = —iX. 



Let us calculate the second derivative of the function j{y) 



sin 9j^\x\ sin 9j^\y\ 
sin 6^ 1 



— 2(|xp + lyp — 1) sin5+|x| sin6'+|?/| + 4|x||?/| cos6'+|x| cos6'+|?/| 
(-cu + z/)3/2(4 + cj-z/)3/2 

8(|x| cos6'+|x| sin6'+|?/| + |y| sin6'+|x| cos6'+|?/|) cos 6*+ ^ 16sin6'+|x| sin6'+|y| cos^6' 



+ 



(-CU + Z/)2(4 + CU-Z/)2 (_^+;y)5/2(4 + ^_^)5/2 

xl cos^^4.|a;| sin6'+|?/| + |?/| sin 6'+|a;| cos6'_|_|?/| sin6'+|a;| sin^^+|y| cos6'4 



(-CJ + Z/)2(4 + CJ-Z/)^ 



(-CJ + Z/)5/2(4 + ,^_^y)5/2 



u; + 2 -z/) 



Since 



then 



and 



sin6'_|_|z| ~ \z\^v — a;, v ^ and sin6'+|2;| ~ |z|-\/4 + a; — z/, ^ a; + 4, 



\r\y)\ < 



C{l + \xm + \y\^) 

{V - Cj)3/2 



ijj < V < ijJ ^ \^ 



\r\y)\ < 



C(l + 1x^(1 + 



u; + 3<z/<u; + 4. 



(4 + a; - z/)3/2 

For y < an identical calculation leads to the same bound. Therefore the operator valued 
function J-'{i') = Fj^(A + 0) — Fj^(A — 0) satisfies the conditions of Lemma 14.11 with a = u, 
6 = 4 + cj, p = 1/2 and B = i3. ■ 
Next we consider the second summand in the RHS of (14. Sp . 

Lemma 4.4. In the situation of Theorem \A.2\ 



j e^\P{X + 0)-P{X-0)) dX = 0{t'''/^), 



(4.6) 



c+uC- 



in the norm B. 
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Proof We consider only the integral over C+ and one component of the matrix P, for 
example, Pu: 

_ (za-2sin^_)e^^+(l^'l+lfl) + (2 sin^_-za)e*''-(l^'l+l2^l) + i/3(e*^-|j'l+*^+l^l -e^^+ls^l+'^-l^l) 
^^^^ '^'^)= 2^a(sin 9+ + sin ^_ ) - 4 sin 9+ sin 9. + - 

Denote ( = —uj — i\, then sin 9+ = y/Cy/^--C, sin 9^ = y/—2u — Ca/4 + 2u + (. The Taylor 
expansion in as C — ^ 0, Im(^ > implies 

PnituJ + <, X, y) = Po + Pi{x, y)C/' + P^ix, y)O(C), 

where \Pj{x,y)\ < Cj{l + + j = 1,2. Therefore, if X = iuj + i( E Cj^- then 

J^{uj + C) = + < + 0) - Pn(zcu + «C - 0) = C(C'/'), C ^ 

in the norm of B. Similarly, differentiating two times the function Pii{X,x,y) in A, we obtain 
that 

r\Lu + c) = -Puituj + 2C + 0) + Pi[{iuj + 2C - 0) = o(C"'/'), C ^ 

in the norm of B. In the same way 

r'(4 + 00-0 = -A"i(^(4 + - C) + 0) + Pi';(^(4 + a; - C) - 0) = 0{C'^'), C ^ 
Therefore, the function ^(z^), = —iX satisfies the conditions of Lemma [4.11 with p = 1/2 and 

B = B. m 

5 Modulation equations 

In this section we present the modulation equations which allow a construction of solutions 
ip{x,t) of equation (ll.ip close at each time t to a soliton i.e. to one of the functions ipuj{,x) in 
the set 5+ US- described in section |2] with time varying ("modulating") parameters {uj,9) = 
{uj{t),9{t)). Let us rewrite f 1 1.1 1) in the real form 

j^{x,t) = -nx,t)-6{x)F{^{0,t)), (5.1) 

as an equation for ip{x,t) G with F{iIj) G which is the real vector version of F{iIj) G C 
Then ^jJ{x, t) = e-'^(*^$^(t)(x) is a solution of (15. ip if and only if ^ = a; and ti; = 0. 
We look for a solution to (15. ip in the form 

ij{x,t) = e^'^W($^(i)(x) + x(x,t)) = e^'^Wv]/(a;,t), ^ii{x,t) = <l>^(i)(x) + x(x, t). (5.2) 

Since this is a solution of (15. ip as long as x = and 9 = uj and ci; = it is natural to look for 
solutions in which x is small and 

9{t) = [ uj{s)ds + -f{t) 



with 7 treated perturbatively. Observe that so far this representation is underdetermined since 
for any (^uj{t),9{t)^ it just amounts to a definition of x; it is made unique by restricting x(t) 
to lie in the image of the projection operator onto the continuous spectrum = P'^{uj{t)) or 
equivalently that 

P°x(t) = 0, P? = P°(^(t)) = / - P'^iioit)) (5.3) 

Now we give a system of modulation equations for uj{t), 7(t) which ensure the conditions (15. 3p 
are preserved by the time evolution. 
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Lemma 5.1. (i) Assume given a solution of Ii5.1\) with regularity as described in Theorem \2.1\ 
which can be written in the form ( (5.2)) -/ (5.3)) with continuously differentiable Lj(t), 9{t). Then 



X = Cx- djd^^^ + 7J"^($.. + X) + Q 



(5.4) 



where Q(x, uj) = -6{x)j~^ (F(<l>^ + x) - F($^) - F'(<l>^)x) , and 



(P°Q,^) 



(5.5) 



(5.6) 



where = P^{uj{t)) is the projection operator defined in / (3.10)) and d^^P^^ = 9^P°(u;) evaluated 
at u = uj{t). 

(a) Conversely given ip a solution of ( f5.1l) as in Theorem \2.1\ and continuously differentiable 
functions uj(t), 6{t) which satisfy ( [5.5j) -( f5l6]) then x defined by ( 15. 2 j) satisfies ( 15.41) and the 
condition ( (5.3|) holds at all times if it holds initially. 

Proof This can be proved as in ^ Prop. 2. 2]. ■ 
It remains to show, for appropriate initial data close to a soliton, that there exist solutions 
to (I5.5p - fl5.6l) . at least locally. To achieve this observe that if the spectral condition {SP} holds 
then by Lemma 12.41 the denominator appearing on the right hand side of (15. 5p and (15.61) does 
not vanish for small ||y||/i. This is because 



as discussed in section |2l This has the consequence that the orthogonality conditions really 
can be satisfied for small x because they are equivalent to a locally well posed set of ordinary 
differential equations for t —>■ {6{t),uj{t)). This implies the following corollary: 

Corollary 5.2. (i) In the situation of (i) in the previous lemma assume that ( (5.7j) holds. If 
WxWf is sufficiently small for some p, P the right hand sides of ( (5.5)) and ( (5.6)) are smooth in 
6,u! and there exists continuous TZ = 7l{uj,x) such that 



(a) Assume given ip, a solution of ( (5.1)) as in Theorem \2.1[ Ifuo satisfies ( (5.7)) and x{x, 0) = 
e~^^"tp{x,0) — ^i^g{x) is small in some K norm and satisfies ( (5.3)) there is a time interval on 



6 Time decay for the transversal dynamics 

Let us represent the initial data ipQ in a convenient form for application of the modulation 
equations: the next Lemma will allow us to assume that ([5.3p holds initially without loss of 
generality. 




(5.7) 



7(t)| <7^|x(0,t)|^ 



u{t)\<n\x{o,t)\\ 
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Lemma 6.1. In the situation of Theorem 12.61 there exists a solitary wave ipcjo satisfying the 
spectral condition {SP} such that in vector form 

and for Xoi^) we have 

P°(cDo)(xo) = 0, (6.1) 

and 

WXoW i^nP = d = 0{d) asd^O. 

Proof This can be proved as in PQ Lemma 10.1] by a standard application of the imphcit 
function theorem. ■ 

In Section [8] we will show that our main Theorem 12.61 can be derived from the following 
time decay of the transversal component x{t)'- 

Theorem 6.2. Let all the assumptions of Theorem \2.6\ hold. Ford sufficiently small there exist 
functions t ^ (^uj(t),'j{t)) defined for t > such that the solution ip{x,t) of lib. I]) can be 
written as in / I5.2H5.3)) with / I5.5H5.6)) satisfied, and there exists a number M > 0, depending 
only on the initial data, such that 

M{T) = sup [(1 + tf'Wxmi-, + (1 + ^)'(l7l + 1^1) I] < M, (6.2) 

0<t<T 

uniformly in T > 0, and M = 0{d) as d —>■ 0. 

Remarks 6.3. (0) This theorem will be deduced from Proposition U .1\ in the next section. 

(i) Theorem U.W implies that the norms in the definition of M are continuous functions of 
time (and so M is also). 

(a) The result holds also for negative time with appropriate changes since ip{x,t) solves 
( II. Ij) if and only if ipi^x, —t) does. 

(Hi) The result implies in particular that t^\9 — c^l + t^\uj\ < C, hence uj(t) and 6{t) — tu!+ 
should converge as t ^ oo while ip{x,t) — e^^^^^^ui{t)ix) have limit zero in /^^(M). 



7 Proof of transversal decay 

7.1 Inductive argument (proof of Theorem 16^2]) 

Let us write the initial data in the form 

M^)=e'''{<^u.oix)+xo{x)). (7.1) 

with d = ||xo||/in//i sufficiently small. By Lemma 16.11 we can assume that P°(c<Jo)(xo) = 
without loss of generality. Then the local existence asserted in Corollary 15.21 implies the 
existence of an interval [0,ti] on which are defined functions t (co'(t), 7(t)) satisfying 
(I5.5l) - fl5.6l) and such that M(ti) = p for some ti > and p > 0. By continuity we can make p 
as small as we like by making d and ti small. The following Proposition is proved in section 
El below. 
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Proposition 7.1. In the situation of Theorem 16.21 let M{ti) < p for some ti > and p > 0. 

Then there exist numbers di and pi, independent ofti, such that 

M{ti) < p/2 (7.2) 

ifd= llxollLin//! < di and p < pi. 

Assuming the truth of Proposition 17.11 for now Theorem 16.21 will follow from the next argu- 
ment: 

Consider the set T of ti > such that (ci;(t), 7(t)) are defined on [0, ti] and M(ti) < p. This set 
is relatively closed by continuity. On the other hand, (17. 2p and Corollary 15.21 with sufficiently 
small p and d imply that this set is also relatively open, and hence supT = +oo, completing 
the proof of Theorem 16. 2[ ■ 

7.2 Frozen linearized equation 

A crucial part of the proof of Proposition 17.11 is the estimation of the first term in M, for 
which purpose it is necessary to make use of the dispersive properties obtained in sections [3] 
andm Rather than study directly (15. 4p . whose linear part is non-autonomous, it is convenient 
(following [21 [3]) to introduce a small modification of (15.21) . which leads to an autonomous 
linearized equation. This new ansatz for the solution is 

ip{x,t) = e^\<l>u;{x) + e-^^^-^'>r]), where e{t) = uit + ^o, = andcji = uj{ti) (7.3) 

so that, T] = e^^^~^^x ^i-nd x = e~^^^^^^r]. Since 

X = e-'^^-^^ (?) - J (^ + 7 - ^i)^) 

equation (15. 4p implies 

f, = j-\cu^ - u;)r] + e^^'-'^C (e-^^^-^'^r/) + e^^'-'^ [r'i^u. - djd^^^ + Q[e-^^'-''^r]]^ . (7.4) 
The matrices C and e^'^, where = ^ — ^, do not commute: 

Ce^^ - e^^C = (5(x)6sin0 a, where = ( J ) ' ^ = (7-5) 
Using (17. 5p we rewrite equation (17.41) as 

= j-\uji - u)7] + Cr] + e^(^"^') (^-6{x)bsm{e - e)ar] + j-^7<l>^ - djd^<^^ + Q[e-^'^^-^^r]]^ . 

To obtain a perturbed autonomous equation we rewrite the first two terms on the RHS by 
freezing the coefficients at t = ti. Note that 

r'(^i - ^) + c = Ci - r's{x){v - Vi), 

where V = a + bPi, V\ = V^(ti), and Ci = C(ti). The equation for rj now reads 

= Ci7]-r^6{x){V-Vi)r]+e^^^-^^(^-6{x)bsm{e-e)ar]+j-^(^-djdu,(^^ + Q[e-^^^^^^ (7.6) 
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The first term is now independent of t; the idea is tliat if there is sufficiently rapid convergence 
of io{t) as t ^ oo the other remaining terms are small uniformly with respect to ti . Finally the 
equation (17. 6p can be written in the following frozen form 

r] = CiV + fi (7.7) 

where 

fi = -j-^6{x){V -Vi)r] + e^^^-^^ (^-6{x)b sm{e - e)a7] + j-^<l> - ijd^<l>^ + Q[e-^^^^^^ (7.8) 

Remark 7.2. The advantage of ( (7.71) over ( (5.41) is that it can be treated as a perturbed au- 
tonomous linear equation, so that the estimates from section^ can be used directly. The addi- 
tional terms in fi can be estimated as small uniformly in ti: see lemma U^ below. This is the 
reason for introduction of the ansatz ( (7.3|) . 

Lemma 7.3. In the situation of Proposition \T7\\ there exists c > 0, independent ofti, such that 
forO<t<ti 

\a{t) - ail + \b{t) - 6i| + \e{t) - e{t)\ < cp, 

where 

p:= sup (l + t3)(|^(t)| + |,0(t)|) <M(ti). (7.9) 

0<t<ti 

Proof By f [7.9p . we have 



\a{t) - a{ti)\ = I / a(r)dr| < c( sup (1 + r^)|u;(r)|) / - 



dr 

— 2 ^ ^Z^' 



since \d{r)\ < c\uj{t)\. The difference \b{t) — b{ti)\ can be estimated similarly. Next 

e{t)-9{t) = [ u;(r)rfr + 7(t) -u;(ti)t- 7(0) = [ {uj{t) - uj{ti))dT + [ ^{T)dT 
Jo Jo Jo 

ft ftl pt 

uj{s)dsdT+ / 7(r)rfr. (7.10) 



Jt 



By f [7.9[) the first summand in RHS of f [7.10[) can be estimated as 

»i rr ft rti 



^ ^ \d;{s)\ds dr < J\l + sf+^\d;{s)\-^^-^^dsdT 



<c sup {1 + sy+'\dj{s)\ [ [ j—^^^dsdT<cp 

0<s<ti Jo Jt (1 + 

since the last integral is bounded for t G [0,ti]. Finally, for the second summand on the RHS 
of O7.10j) inequality f[7.9j) imphes 

I / 7(r)rfr|<c sup {1 + T^)\'y{T)\ / — — ^ < cp 
Jo o<T<ti J 1 + 
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7.3 Projection onto discrete and continuous spectral spaces 

From sections [3] and m we have information concerning U{t) = e^^*, in particular decay on the 
subspace orthogonal to the (two dimensional) generalized null space. It is therefore necessary 
to introduce a further decomposition to take advantage of this. Recall, by comparing fl5.2p and 
fra that 

r] = e^('-^^X and P?x(t) = 0. (7.11) 

Introduce the symplectic projections = P° and P^ = Pj^ onto the discrete and continuous 
spectral subspaces defined by the operator Ci and write, at each time t G [0,ti]: 

7^{t) = g{t) + h{t) (7.12) 

with g{t) = P5?7(t) and h{t) = Plr]{t). The following lemma shows that it is only necessary to 
estimate h(t). 

Lemma 7.4. In the situation of Proposition ITTTl assume 

sup (|c^(t) - cuil + \9{t) - 9i{t)\) = A 

0<t<ti 

is sufficiently small. Then for <t <ti there exists c(A,co'i) such that 

c(A,c^i)~^||/i||;oo^n/2 < Ihlli^^nP < c(A,c^i)||/i||;oo^n/2. (7.13) 
Proof This can be proved as in ^ Lemma 11.5] ■ 

7.4 Proof of Proposition 17.11 

To prove Proposition 17.11 we explain how to estimate both terms in M, (16. 2p . to be < p/4, 
uniformly in ti. 

Estimation of the second term in M. As in Corollary 15.21 we have 

M(t)'^ 

i7(t)i + \m\ < coix(o,t)r < ^0(1:^^' ^ ^ 

since |x(0,t)| < ||x(^)lbf° • Finally let pi < l/(4co) to complete the estimate for the second 
term in M as < p/4. 

Estimation of the first term in M. By Lemma [7.41 it is enough to estimate h. Let us apply the 
projection P^ to both sides of (17. 7p . Then the equation for h reads 

/i = Ci/i + P^fi (7.14) 

Now to estimate h we use the Duhamel representation: 



h{t) = U{t)h{0)+ [ U{t - s)PlU{s)ds, t<ti. 
Jo 



(7.15) 



with U(t) = e*^^* the one parameter group just introduced. Recall that Plh(t) = for t G [0, ti]. 
Therefore 

||f/(t)MO)||,^^ <c(l + t)-^/^||MO)||,.nP <c(l + t)-^/^h(0)||,i^ (7.16) 
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by Theorem 14.21 and inequalities fl3.2p and (17.131) . Let us estimate the integrand on the right- 
hand side of (17.151) . We use the representation (17.81) for fi and apply Theorem 14.21 Corollary 
15.21 and Lemma 17.31 to obtain that 

\\Uit-s)PtU\\i^^ < c(l + t-.)-^/2||P^(fi(t))||,i 

< cii+t-sr'/'(^mt)\'+Mo,t)\ 

< c(l + t-s)-3/2(^||r^(t)||2^^^+p||7;(t)||,^J, t<h. (7.17) 
Now (17T31) . (17T5I) . (061) and (171711 imply 

t 

\\vmiT,<c{i+tr'/^vmiini^+ci f jrT^^-x^.lwvmh+pMsnr: 

Multiply by (1 + tf/^ to deduce 







l+t-s)3/2 V'"'^'"-^ r,r,v'^/,K_, 



t 





t 



I (l+t)3/2(l + s)-3/2 
+ / (l+t-.)3/2 (1 + ^) ^ 11^(^)11'-/^ 



since ||?7(0)||;in;2 < Introduce the majorant 



m{t) :=sup(l + s)=^/2 11^(^)11 ^<^^ 

[0,i] 

and hence 

mit) <cd + c,m\t) j ' ds + pc,m(t) j ^'^,'^^7-^)31''' ^^'^^^ 



It easy to see (by sphtting up the integrals into s < t/2 and s > t/2) that both these integrals 
are bounded independent of t. Thus (I7.19P implies that there exist c, C2,C3, independent of ti, 
such that 

m{t) < cd + pc2m{t) + c^m^{t), t < ti. 

Recall that m{ti) < p < pi by assumption. Therefore this inequality implies that m(t) is 
bounded for t < ti, and moreover, 

m(t) < c^d, t <ti 

if d and p are sufficiently small. The constant C4 does not depend on ti. We choose d in (12.61) 
small enough that d < p/{4c4). Therefore, 

[0,ti] 

if d and p are sufficiently small. This bounds the first term as < p/4 by (17.111) and hence 
M{ti) < p/2, completing the proof of Proposition 17. 1[ ■ 
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8 Soliton asymptotics 

Here we prove our main Theorem 12.61 using the bounds ( 16. 2p . For the solution ilj{x,t) to ( 11. ip 
let us define the accompanying soliton as s{x,t) = 4'oj(t)ix)e'^^^^\ where 9{t) = uj{t) + 'y{t). Then 
for the difference z{x,t) = 'ip{x,t) — s{x,t) we obtain easily from equations ( 11.11) and ( 11.61) 



iz{x, t) = -z"{x, t) + 'ys{x, t) - iud^s{x, t) - 6{x) yF{ip{0, t)) - F(s(0, t))j . (8.1) 

Then 

t 

z{t) = W{t)z{0) + jw{t- r)f{; r)dT, (8.2) 



where f{x, t) = 'js{x, t) —iujd^s{x, t)—6{x) (^F{ip{0, t)) — F(s(0, t)) j , and W{t) is the dynamical 

group of the free Schrodinger equation. Since 7(t) — 7+, uit) — u!+ = 0(t^^), and therefore 
9{t) — u+t — 7+ = 0{t^^) for t — >• 00, to establish the asymptotic behavior (12. 7p it suffices to 
prove that 

z(t) = iy(t)$+ + r+(t) (8.3) 

with some $+ G P{Z) and ||r+(t)|| = 0{t-^/^). Let us rewrite (ESD as 

00 t 
z{t) = W{t) (z(0) + j W{-T)f{; r)dr) -Jw{t- r)f{-, r)dr = (t)$+ + r+(t). (8.4) 


Let us recall that by (16. 2p 

\m\<c{l + tY\ |7(t)| <c(l+t)-3, |F(V^(0,t))-F(s(0,t))| <c|x(0,t)| <c(l+t)-3/2. 

00 

Hence, the unitarity in of the group W{t) implies that = 2;(0) + j W {—t) f {■ ^ T)dT G 



and ||r+(t)|| = 0{t-^l^), t ^ 00. ■ 

A The resolvent 

A.l Calculation of the matrix kernel 

Here we will construct matrix kernel of the resolvent R(A) explicitly 

R(A,x,,)=f I'f.^'^'y] (a.i) 

^ ' '^^ V ^2l(A,X,y) i?22(A,X,y) ) ^ ' 

which is the solution to the equation 

{C-X)K{\x,y) = 5{x-y)(\ \ Y (A. 2) 
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Calculation of first column For the first column Rj{X,x,y) :- 
matrix R(A, x, y) we obtain 



Ru{X,x,y) 
R2ii\x,y) 



{C-X)Rj{X,x,y) = 6{x-y) 



of tlie 



(A. 3) 



If X 7^ and x ^ y, (lA. 3p takes tlie form 



-A 
Di 



D, 



-A 



Ri{\x,y) = I ^ 



-A -Ar-CJ 



-A 



Ri{\,x,y) = 0, xj^O, x^y. 

(A. 4) 

Tlie general solution is a linear combination of exponential solutions of type e^^^v. Substituting 
into (IA7^ . we get 

-A 2-2cos^ + cu \ n 
„ . ]v = 0. (A. 5) 

-2 + 2 cosfc' — "-^ / 

For nonzero vectors v, the determinant of the matrix vanishes, 

A2 + (2-2cos0 + cj)2 = O. 
Finally, we obtain four roots ±^±(A) in D := — tt < Re^^ < tt with 

2 - 2cos^±(A) = -cu T^A. (A. 6) 



We choose the cuts in the complex plane A: 

the cut C4. := [iu, i{uj + 4)] for 6'+(A) and the cut C_ := [— ^(ci; + 4), —iuj] for ^^-(A) and 

lm^±(A)>0, AgC\C±. 
It remains to derive the vector v = [vi, V2) which is solution to flA. SD : 

2-2cos6'± + cj ±iiA 



(A. 7) 



V2 



A 



-Vi = —;—Vi = ±iVi. 



Therefore, we have two corresponding vectors v± 



A 



±1 



and we get four linearly indepen- 



dent exponential solutions f+e^*^+^ and w_e^*^~^. 

Now we can solve the equation (]A. 3\} . First we rewrite it using the representation (13. ip for 
the operator C, 



-A -Al + cjW Rn{X,x,y) 
Al-uj -X J \ R2i{X,x,y) 



6{x-y) 



-6{x) 



a\ / i?n(A,0,y) 
-a-b Oj\ i?2i(A,0,y) 

(A. 



Let us consider y > for the concreteness. Then the RHS vanishes in the open intervals 
(—00, 0),(0, y) and {y, 00). Hence, for the parameter A outside the cuts C±, the solution admits 
the representation 



Ri{X,x,y) = < 



^-g-ie+a-^^ + BZe-^^-^'v- + 5+e*^+^t;+ + Bte'^-^'v-, < x < y, (A. 9) 



V- 



x> y 



16 



since by (lA. 7\} . the exponent e~*^*^ decays for x — oo, and similarly, e'^*^ decays for x 
oo. Next we need eight equations to calculate the eight constants . . . We have two 
continuity equations and two jump conditions for the derivatives at the points x = and x = y. 
These four vector equations give just eight scalar equations for the calculation. 
Continuity at x = y: Ri{y — 0,y) = Ri{y + 0,y), i.e. 

BZv^/e^ + BZv-je^ + i?^f_|_e+ + B^v^e^ = C+v+e^ + C_f_e_, 

where e± := 6*^="=^. It is equivalent to 

5;/e+ + 5+e+ = C+e+, 

5Z/e_ + Bte- = C_e_. 
Continuity at x = 0: Ri{—0,y) = Ri{+0^y), i.e. 

+ A^v^ = i?7?;+ + BZv- + B^v^ + B^v. 



that is equivalent to 



"Jump" at X = y: Noting that 



A_ = BZ+ Bt. 



^ 

Kt X = y equation (lA. Sp reads 

-AC+e+ - iB^e'^+/e+ - iB+e-'^+e+ - iC+e+e'^+ + (2 + uj)iC+e+ = 1/2, 

-AC_e_ + iBZe'^-/e^ + zEle-^'^e^ + zC_e_e^^- - (2 + cu)iC_e_ = 1/2. 
Also we need equation (]A. 8|) at x = y + 1: 



-AC_e*^-(^+i) + iC_e^^-^ + zC_e^^-(^+2) - (2 + cj)iC_e^^-(?^+i) = 0. 
Hence 

-\C+ - iC+e'^+ + (2 + cu)zC+ = iC+e-'^+, 

-AC_ + iC_e*^- - (2 + uj)iC_ = -iC.e-'^-. 
Combining flA. 13p and flA. 14p . we get 

-i5;e^^+/e+ - i5+e-*^+e+ + iC+e+e^-'^+ = 1/2, 

iBZe'^~ /e_+iBte-'^e^ - iC^e.e^'^- = 1/2. 
After substituting of C± from (]A. lOp . the constants B^ cancel and we get 



5; — (-e'^+ + e-'^+) = 1/2, BZ — ie''^- - e-'^~) = 1/2 



e_i_ e_ 



(A. 10) 



(A. 11) 



(A. 12) 



(A. 13) 



(A. 14) 
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and then 



4 sin ^4-' 



B 



4 sin 9^ 



"Jump" at X = 0: Noting that 

a 
-a-b 



a 
-a-b 



-a — b 
-a — b 



• ^ • / ^ • ^ • 



b ^ b 
where a = a + p = -. Hence, 



a 
-a-b 



At a; = equation (1A. 81) reads 

- iA+e'^+ - iB^e-'^+ - iB+e'^+ + (2 + uj)iA+ = i{A+a + 

-XA_ + iA.e'^- + iBZe-'^+ + iB^e'^- - (2 + uj)iA_ = -i{A+(3 + A_a) 
Also we need equation flA. 8\i at x = y -\- 1: 

-\A+e'^+ - iA+ - iA+e^'^+ + (2 + uj)iA+e'^+ = 0, 

-XA^e'^- + iA_ + iA.e^^^- - (2 + uj)iA_e'^- = 0. 

Therefore 



-\A+ - iA+e"^^+ - iA+e''^+ + (2 + uj)iA^ 



0, 



0. 



-AA_ + iA_e-^^- + iA_e^^- - (2 + uj)iA 
Substituting flA. 17p into flA. 16p we get 

-iB^e~'^+ - iBXe'^+ + iA+e-'^+ = i{A+a + A_p) 

iBZe-'^+ + iBte'^- - iA^e~'^- = + A_a) 

Substituting here ( lA. lip , we get after cancellations, 

(e*^+ - e-'^+ + a)BX + ^Bt = -aB^ - (3BZ 

13 BX + (e^^- - e-'^- +a)Bl = -/SB' - aBZ 
Hence, the solution is given by 

Bi 



b: 



1 / 2z sin 9^ + a -(3 
15 \ -(3 2zsin6'+ + a 



a (3 
(3 a 



B: 



where D is the determinant 



D := {2i sin 9+ + a) {2i sin 9. + a) - 
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and , B_ are given by (]A. 15|) . The formulas (]A. 151) and (lA. 18|) imply 
1 / 2ia sin 6'_ + - 



5 



2L) 



2sin^. 



■o \ n+ 1 2ia sin ^+ + _ 
e+ + z/?e_), 5+ = — f-z/?e+ + ^^-^ 



2 sine. 



Using the identities 



(A. 20) 



2ia sin e_ + — = D — 2ia sin 6*+ + 4 sin 9+ sin 6'_ , 
2ia sin e+ + — = D — 2ia sin 6- + A sin 6*+ sin 6'_ , 



let us rewrite ( lA. 201) as 



5J 



4sine+ 2D 



^(^(m-2sine_)e++z/5e_), 5^ 



1 



4sine_ 2D 



i[3e^+{ia—2 sin6'+)e_ 



(A. 21) 

Finally, the formulas (lA. 9I) - (]A. 11|) . (]A. 15^ and (lA. 2ip give the first column Rj{\,x,y) of 
the resolvent for y > 0: 

RiiX, X, y) = r,(A, X, 2/) + P,(A, 2/), (A. 22) 

where 



r^fA, X, y) = (e'^+\^-y\ - e*''+(l^'l+IJ^l))t;+ (e'^-\-y\ - e*«-(l-l+ls/l))^;„, (A. 23) 

4 sin 64. 4 sin 6^ 



and 



Pi{\x,y) 



2D L 



(ia - 2sine_)e*^+(l^'l+l^l) + i/5e*(^+l^l+^-l^l))t;+ (A. 24) 



i/3e 



i(9.\x\+9+\y\) 



+ (m-2sine+)e*^-(l^l+l^l) 



Calculation of second column The second column is given by similar formulas with the 
vector I ^ j instead of f ^ j in (]A. Sp . Respectively, (lA. 12^ is changed by 



Hence, we have now change 1/2 by —i/2 in the first equation of (]A. lOP and 1/2 by i/2 in the 
second one. Respectively, ( lA. 151) for the second column reads 



B- 



B- 



4 sin e+ 4 sin 6*. 

Then the second column Rij{\,x,y) of the resolvent reads: 

Rji{X,x,y) = Tjj{X,x,y) + Pji{X,x,y), 

where 



Tii{X,x,y) 



4 sin 9 



! (^iO+lx-yl _ ^i9+i\x\ + \y\)^^^ _ ^ (^id.\x-y\ _ ^ie ^ {\x\ + \y\) 



and 



Pii{X,x,y) = — 



4 sin 9_ 

ia - 2sine_)e^^+(l^l+l^l) + i/^e^^^+l^l+^-l^l))^;^ 
+ [^f3e'^s.\x\+eM) _ - 2sine+)e^^-(l"l+l^l))t; 



(A. 25) 



)v-, (A. 26) 



(A. 27) 



Note, that if ?/ < we get the same formulas. 
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A. 2 The poles of the resolvent 

The poles of the resolvent correspond to the roots of the determinant (lA. 19^ . 

D{\) ■.= a'^ + 2m(sin 9+ + sin O.) - 4 sin 9+ sin 9^ - = 0. (A. 28) 

with 9± as in (]A. 6|) - (1A. 7\i . Thus -D(A) is an analytic function on C \ C_ U C+. Since there are 
two possible values for the square roots in 9± there is a corresponding four-sheeted function 
D{\) analytic on a four sheeted cover of C which is branched over C_ and C_|_. We call the sheet 
defined by (1A. 71) the physical sheet. 

We will reduce the equation (lA. 281) to the solution of two successive quadratic equations. 
These can be solved explicitly but the process involves squaring and thus actually produces 
zeros of the function -D(A) rather than of -D(A). Therefore we will then have to check whether 
or not the roots do actually lie on the physical sheet. 

Proposition A.l. If a' = (4a + a^)/4C^ then X = is a root of the determinant D{X) with 
multiplicity 4, otherwise X = is a root of the determinant -D(A) with multiplicity 2. 

Proof First let us check that A = is a root of -D(A). Let us represent uj by means of 
a = a{C^). By resuhs of g2] 

cosh k = — — — ^ , sinh /c = — . 

2 ' 2 

Therefore 

(cj + 2)2 = 4 + a^ (A. 29) 

D{0) =a'^-p^ + 2iaia + = (a + h/2f - I? /A - 2(a + h/2)a + = 

since 

sin^i = - + 2)2/4 = - (4 + a2)/4 = ^-a2/4 = ia/2. 

Now let us compute -D'(A): 

. / -i{uj + 2 + iX) i{uj + 2-iX) 



+ 



V4 - (cu + 2 + 2A)2 ^/4^r(JT2^^iW' 
i{uj + 2 + iX)>yA -{uj + 2- iXy i{uj + 2- iX)>yA - {u + 2 + iXy 
^A-{uj + 2 + iA)2 ^4 - (cj + 2 - iXy 



Hence -D'(O) = and A = is the root of -D(A) of multiplicity at least 2. Further calculation 
shows that the Taylor series for D near zero takes the form: 



D(A) = + -^\ ^ + 0{X% (A. 30) 



Therefore A = is the root of -D(A) of multiplicity 4 if and only if 6 = (4a + a^)/2, i.e. 
a' = (4a + a3)/4C2. ■ 
Now we prove that there exist a values of a'C^ = f{a) such that the only root of D{X) is 
A = with multiplicity 2. 

Proposition A. 2. If a' = —a/C^, then the only root of D{X) is X = with multiplicity 2. 
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Proof One has b = —2a, a = a + | = and then equation (]A. 281) imphes 

sin e+ sin 9_ = -p'^/A = -a^/A < 0. (A. 31) 

From the other hand, by flA. 29^ we get 

sm.^ 9+ sm^ 9_ = yl - —j - —j (A. 32) 

2 \2 



\^-- — T^^-r) ^- 7^ — = + ^ +^ +^ 

V 4 4/ 4 V44/ 4 



Hence, (lA~3T|l and (1X7121) imply 

^4 , „2 \2. 2 ^2 \2 

I6 = (-t + t) +^ + — 

Hence 

A2(A2 + 2a2 + 16) = 



It remains to check that the roots A = ■^i\/2a^ + 16 don't he on the physical branch. 
Let 6^+ = x+ + 2?/+, = a;_ + iy^ with ?/± > 0. Then 

^ + 2 + V2a2 + 16 aAM^ + ^2a2 + 16 ^ 
cos = COS X4- cosh y+ — 1 sm x+ smh ?/+ = ^ = ^ > U, 

oj + 2- yj2a? + 16 \/a2T4 - V2a2 + 16 ^ 

cos y_ = cos x_ cosh y_ — I sm a;_ smh y_ = = < 0. 

^ ^ 2 2 

Hence, x+ = 0, x_ = vr and 

sin 6'+ = sin(iy+) = isinh?/+, sin^_ = sin(7r + iyS) = — isinh?/_ 
Finally we obtain that sin6'+sin6'_ = sinh sinh ?/„ > 0, which contradict flA. 3ip . I 



References 

[1] V.Buslaev, A.Komech, E.Kopylova, D.Stuart, On asymptotic stability of solitary waves in 
Schrodinger equation coupled to nonlinear oscillator. Comm. Partial Diff. Equations 33 
(2008), 669-705. 

[2] V.S. Buslaev, G.S. Perelman, On the stability of solitary waves for nonlinear Schrodinger 
equations, j4mer. Math. Soc. Trans. (2) 164 (1995), 75-98. 

[3] V.S. Buslaev, C. Sulem, On asymptotic stability of solitary waves for nonlinear Schrodinger 
equations, 74nn. Inst. Henri Poincare, Anal. Non Lineaire 20(2003), no. 3, 419-475. 

[4] P.A. Deift, A.R. Its, X. Zhou, Long-time asymptotics for integrable nonlinear wave equa- 
tions, pp. 181-204 in: A.S. Fokas,V.E. Zakharov, (ed.). Important Developments in Soliton 
Theory, Springer, Berlin, 1993. 

[5] M. Grillakis, J. Shatah, W.A. Strauss, Stability theory of solitary waves in the presence of 
symmetry, I; II. J. Func. Anal. 74 (1987), no.l, 160-197; 94 (1990), no.2, 308-348. 



21 



[6] V. Imaikin, A.I. Komech, B. Vainberg, On scattering of solitons for the Klein-Gordon 
equation coupled to a particle, submitted to Comm. Math. Phys., 2005. 

[7] A. Jensen, T. Kato, Spectral properties of Schro dinger operators and time-decay of the 
wave imiciioYis^Duke Math. J.46 (1979), 583-611. -307. 

[8] A. Komech, E. Kopylova, Scattering of solitons for Schrodinger equation coupled to a 
particle, Russian J. Math. Phys. 13 (2006), no. 2, 158-187. 

[9] E.A. Kopylova, On existence of solitary waves in discrete Schro dinger equation coupled to 
a nonlinear oscillator, submited to Russian J. Math. Phys. 

[10] J. Miller, M. Weinstein, Asymptotic stability of sohtary waves for the regularized long-wave 
equation. Comm. Pure Appl. Math. 49 (1996), no. 4, 399-441. 

[11] R.L. Pego, M.I. Weinstein, On asymptotic stability of solitary waves, Phys. Lett. A 162 
(1992), 263-268. 

[12] R.L. Pego, M.I. Weinstein, Asymptotic stability of solitary waves, Commun. Math. Phys. 
164 (1994), 305-349. 

[13] C.A. Pillet, C.E. Wayne, Invariant manifolds for a class of dispersive, Hamiltonian, partial 
differential equations, J. Differ. Equations 141 (1997), No. 2, 310-326. 

[14] A. Soffcr, M.I. Weinstein, Multichannel nonlinear scattering for nonintegrable equations. 
Comm. Math. Phys. 133 (1990), 119-146. 

[15] A. Soffer, M.I. Weinstein, Multichannel nonlinear scattering for nonintegrable equations. 
II. The case of anisotropic potentials and data, J. Differential Equations 98 (1992), no. 2, 
376-390. 

[16] A. Soffer, M.I. Weinstein, Selection of the ground states for NLS equations. Rev. Math. 
Phys. 16 (2004), no. 8, 977-1071. 

[17] D.M.A. Stuart, Modulational approach to stabihty of non-topological solitons. Journal de 
Mathematiques Pures et Appliqu'ees 80 (2001) , no. 1, 51-83. 

[18] M. Weinstein, Modulational stabihty of ground states of nonlinear Schrodinger equations, 
SI AM J. Math. Anal. 16 (1985), no. 3, 472-491. 

[19] A. Zygmund, Trigonometric Series I, Cambridge, 1968. 



22 



